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Abstract 

Three decades ago Heath found the integral form of the exact an- 
alytic growing mode solution of linear density perturbation 5 in sub- 
horizon scales including the cosmological constant or the curvature 
term. Interestingly, we are able to obtain the analytic solution for gen- 
eral dark energy models with the constant equation of state u)de- We 
compare the correct analytic growing mode solution S with the text 
book solution S D . Indeed, both solutions are equal to each other when 
bJde = — 1 • We also able to extend this solution for the specific form of 
time varying u)de- 

The background evolution equations in a fiat Friedmann-Robertson- 
Walker universe (p m + pd e = Per) are 
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where uide is the equation of state (eos) of dark energy, p m and pd e are the 
energy densities of the matter and the dark energy, respectively. We consider 
the constant oj^e- The sub-horizon scale linear perturbation equation with 
respect to the scale factor a are given in the reference [1], 
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We rewrite the above equation 
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where x = In a and = (^-a 3tJde + 1) 1 = (Y + We are able to find 
the exact analytic growing mode solution of 5 for any value of the constant 
ujde- After replacing new parameter Y in the equation (4), we have 
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Now we try S(Y) = cY a B(Y) because it is the most general combination of 
the solution for the above equation (5). Now we replace 5 into the above 
equation (5) to get, 
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There are two alternative ways to make the above equation as the Hyperge- 
ometric differential equation, Y = —X or 1 + Y = X. The complete solution 
of the above equation becomes 
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where F is the hypergeometric function. Thus, the full analytic solution of 
the sub-horizon scale linear perturbation becomes 
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This analytic solution does not have any physical meaning before we fix 
the coefficients c\ and C2- If we want to have the correct growing mode 
solution from the above analytic solution, then this solution should follow 
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Figure 1: Evolutions of 5 and 5 D . a) When oj^e = —0.9. b) For u^ e = — 1. 



the behavior of growing mode solution at early epoch a ~ 0.1. In other 
world, the coefficients of the 5 should be fixed by using the initial conditions 
for the growing mode solution 
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We compare this correct analytic solution (5) with the textbook solution 
(S D ) given in the reference [2] when we use the definition of growing mode 
as in the text book [3]. 
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In Fig. la, we compare the behavior of them. The horizontal axis is scale 
factor a and the vertical line indicated the amplitude of 5. The solid line 
is the correct growing mode solution 5 and the dashed line is for 5 D . In 
this figure we use the initial conditions (9) to get the value of c\ and 02. 
For example, (c\,C2) = (1.088, —1.223) for 5. Definitely, 5 is bigger than S D 
because S D was suppressed by the additional source term in the perturbation 
equation [2]. In Fig. lb, we show the behaviors of 8 and S D when code = — 1. 
Definitely, two solutions are exactly matched to each other, even though 
the formula for two solutions look quite different for co^e = — 1- From the 
initial conditions we find (01,02) = (cf,^) = (1.085,-0.943), which are 
used for the figure lb. We can also find the behavior of the decay mode 
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Figure 2: a) Evolution of 5 for different values of ujde- b) Evolutions of /. 
solution 5d from this analytic solution after we choose the c\ and C2 by using 

—3/2 

the another initial conditions for decay mode solution 5d(ai) oc a i and 

Itola, oc — fa, 5 ^ 2 . Thus, we can find the both growing and decaying mode 
solutions from this analytic solution 5 without any ambiguity. 

We need to check the growth index and growth index parameter based 
on this correct linear perturbation equation. Definitely, the values of these 
quantities are changed if we use the correct growing mode solution [4]. 5 D 
has the additional source term —AttG(1 + ujde)(l + ^de)Pde- $ D oscillates 
at late time when — 1 < code < — 1/3 by including this term. Because this 
additional source term becomes positive and acts as a restoring force. Thus, 
the growth index or growth index parameters obtained from this solution 
are bigger than the correct values. For u>de < — 1 or u>de > —1/3, this 
source term has the negative sign and give the additional contribution to 
the correct source term which should be just from the matter —4irGp m if we 
assume that the dark energy is homogeneous. Thus, in this case 5 D grows 
faster than the correct solution 5 and gives the smaller value of growth index 
parameter than the correct one. 



We show the behavior of 5 for the different values of Ude in Fig. 2a. Dotted, 
solid, and dashed lines correspond to Ude = —1.2,-1.0 and —0.8, respec- 
tively. We have more matter ratio in the past for the smaller values of ujde 
to give the larger values of 5 at present. The evolutions of the growth index 
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are depicted in Fig. 2b with the same notations as Fig. 2a. 
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Figure 3: The evolution of code m Eq. 11. 



We are able to extend this analytic solution in the specific form of time 
varying uj^e ■ We choose the parametrization of o^e 
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. We show the behavior of this parametrization in one specific case, u)\ = 
and u)2 = 0.01 in Fig. 3. It rapidly changes from z = 1 and reaches to 
— 1 at near present. Definitely, it diverges at present. One can change the 
slope of evolution, the present and past values from the proper values of U\ 
and L02- We want to show the possibility of extension of the exact analytic 
solution of the growth factor to the time varying ujd e - We will not deep into 
the detail or the validity of this parametrization at this moment. 



With the parametrization of Eq. (11), we are able to find the exact analytic 
solution of S. 
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where X = Qe 3uJ2 a 3uJ1 and Q = . It is straight forward to get this solution 

de 

from the cvde parametrization given n Eq. (11). 
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The exact analytic solution provides the convenient and economic tools 
for probing the properties of sub-horizon scales growth factor and observa- 
tional quantities related to it [4]. 

We thanks Y. Gong, E. Linder A. Starobinsky, S. Varun, and A. Wang 
for useful comment, K.-W. Ng for fruitful discussion. We especially thanks 
S. Habib for the correct reference and pointing out the mistake in the text- 
book solution. 
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